Genetic variant of the sheep E2F8 gene and its associations with litter size.

The economic efficiency of sheep breeding, aiming to enhance productivity, is a focal point for improvement of sheep breeding. Recent studies highlight the involvement of the Early Region 2 Binding Factor transcription factor 8 (E2F8) gene in female reproduction. Our group's recent genome-wide association study (GWAS) emphasizes the potential impact of the E2F8 gene on prolificacy traits in Australian White sheep (AUW). Herein, the purpose of this study was to assess the correlation of the E2F8 gene with litter size in AUW sheep breed. This work encompassed 659 AUW sheep, subject to genotyping through PCR-based genotyping technology. Furthermore, the results of PCR-based genotyping showed significant associations between the P1-del-32bp bp InDel and the fourth and fifth parities litter size in AUW sheep; the litter size of those with genotype ID were superior compared to those with DD and II genotypes. Thus, these results indicate that the P1-del-32bp InDel within the E2F8 gene can be useful in marker-assisted selection (MAS) in sheep.

Deep whole-genome analysis of 494 hepatocellular carcinomas.

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

Whole-Genome Sequencing Reveals Rare Off-Target Mutations in MC1R-Edited Pigs Generated by Using CRISPR-Cas9 and Somatic Cell Nuclear Transfer.

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been widely used to create animal models for biomedical and agricultural use owing to its low cost and easy handling. However, the occurrence of erroneous cleavage (off-targeting) may raise certain concerns for the practical application of the CRISPR-Cas9 system. In this study, we created a melanocortin 1 receptor (MC1R)-edited pig model through somatic cell nuclear transfer (SCNT) by using porcine kidney cells modified by the CRISPR-Cas9 system. We then carried out whole-genome sequencing of two MC1R-edited pigs and two cloned wild-type siblings, together with the donor cells, to assess the genome-wide presence of single-nucleotide variants and small insertions and deletions (indels) and found only one candidate off-target indel in both MC1R-edited pigs. In summary, our study indicates that the minimal off-targeting effect induced by CRISPR-Cas9 may not be a major concern in gene-edited pigs created by SCNT.

Systematic Comparison of Computational Tools for Sanger Sequencing-Based Genome Editing Analysis.

Successful genome editing depends on the cleavage efficiency of programmable nucleases (PNs) such as the CRISPR-Cas system. Various methods have been developed to assess the efficiency of PNs, most of which estimate the occurrence of indels caused by PN-induced double-strand breaks. In these methods, PN genomic target sites are amplified through PCR, and the resulting PCR products are subsequently analyzed using Sanger sequencing, high-throughput sequencing, or mismatch detection assays. Among these methods, Sanger sequencing of PCR products followed by indel analysis using online web tools has gained popularity due to its user-friendly nature. This approach estimates indel frequencies by computationally analyzing sequencing trace data. However, the accuracy of these computational tools remains uncertain. In this study, we compared the performance of four web tools, TIDE, ICE, DECODR, and SeqScreener, using artificial sequencing templates with predetermined indels. Our results demonstrated that these tools were able to estimate indel frequency with acceptable accuracy when the indels were simple and contained only a few base changes. However, the estimated values became more variable among the tools when the sequencing templates contained more complex indels or knock-in sequences. Moreover, although these tools effectively estimated the net indel sizes, their capability to deconvolute indel sequences exhibited variability with certain limitations. These findings underscore the importance of judiciously selecting and using an appropriate tool with caution, depending on the type of genome editing being performed.

Relationship between Angiotensin-Converting Enzyme Insertion/Deletion Polymorphism and the Risk of COVID-19: A Meta-Analysis.

Introduction: Research shows the correlation between angiotensin-converting enzyme (ACE) deletion and insertion (D/I) polymorphism and COVID-19 risk; yet, conclusive evidence is still lacking. Thus, a meta-analysis of relevant articles was performed to more accurately estimate the relationship of ACE I/D polymorphism with the risk of COVID-19. Material and Methods. Relevant literature from the PubMed database was systematically reviewed, and odds ratios (ORs) and associated 95% confidence intervals (CIs) were measured. Additionally, the metapackage from Stata version 15.0 was used for statistical analysis. Results: The meta-analysis eventually contained 8 studies, including 1362 COVID-19 cases and 4312 controls. Based on the data, the ACE I/D polymorphism did not show an association with COVID-19 risk (D vs. I: OR = 1.25, 95% CI = 0.96-1.64; DD vs. II: OR = 1.89, 95% CI = 0.95-3.74; DI vs. II: OR = 1.75, 95% CI = 0.92-3.31; dominant model: OR = 1.88, 95% CI = 0.99-3.53; and recessive model: OR = 1.24, 95% CI = 0.81-1.90). Further, subgroup analyses stratified based on case proved that the ACE D allele demonstrated an association with increasing risk of COVID-19 severity (D vs. I: OR = 1.64, 95% CI = 1.01-2.66; DD vs. II: OR = 4.62, 95% CI = 2.57-8.30; DI vs. II: OR = 3.07, 95% CI = 1.75-5.38; dominant model: OR = 3.74, 95% CI = 2.15-6.50; and recessive model: OR = 1.28, 95% CI = 0.46-3.51). Conclusions: The ACE D allele was clearly related to an enhanced risk of COVID-19 severity. Hence, it is imperative to take into account the influence of genetic factors during the development of future vaccines.

The landscape of genomic structural variation in Indigenous Australians.

Indigenous Australians harbour rich and unique genomic diversity. However, Aboriginal and Torres Strait Islander ancestries are historically under-represented in genomics research and almost completely missing from reference datasets1-3. Addressing this representation gap is critical, both to advance our understanding of global human genomic diversity and as a prerequisite for ensuring equitable outcomes in genomic medicine. Here we apply population-scale whole-genome long-read sequencing4 to profile genomic structural variation across four remote Indigenous communities. We uncover an abundance of large insertion-deletion variants (20-49 bp; n = 136,797), structural variants (50 b-50 kb; n = 159,912) and regions of variable copy number (>50 kb; n = 156). The majority of variants are composed of tandem repeat or interspersed mobile element sequences (up to 90%) and have not been previously annotated (up to 62%). A large fraction of structural variants appear to be exclusive to Indigenous Australians (12% lower-bound estimate) and most of these are found in only a single community, underscoring the need for broad and deep sampling to achieve a comprehensive catalogue of genomic structural variation across the Australian continent. Finally, we explore short tandem repeats throughout the genome to characterize allelic diversity at 50 known disease loci5, uncover hundreds of novel repeat expansion sites within protein-coding genes, and identify unique patterns of diversity and constraint among short tandem repeat sequences. Our study sheds new light on the dimensions and dynamics of genomic structural variation within and beyond Australia.

VarSCAT: A computational tool for sequence context annotations of genomic variants.

The sequence contexts of genomic variants play important roles in understanding biological significances of variants and potential sequencing related variant calling issues. However, methods for assessing the diverse sequence contexts of genomic variants such as tandem repeats and unambiguous annotations have been limited. Herein, we describe the Variant Sequence Context Annotation Tool (VarSCAT) for annotating the sequence contexts of genomic variants, including breakpoint ambiguities, flanking bases of variants, wildtype/mutated DNA sequences, variant nomenclatures, distances between adjacent variants, tandem repeat regions, and custom annotation with user customizable options. Our analyses demonstrate that VarSCAT is more versatile and customizable than the currently available methods or strategies for annotating variants in short tandem repeat (STR) regions or insertions and deletions (indels) with breakpoint ambiguity. Variant sequence context annotations of high-confidence human variant sets with VarSCAT revealed that more than 75% of all human individual germline and clinically relevant indels have breakpoint ambiguities. Moreover, we illustrate that more than 80% of human individual germline small variants in STR regions are indels and that the sizes of these indels correlated with STR motif sizes. VarSCAT is available from https://github.com/elolab/VarSCAT.

A CRISPR-del-based pipeline for complete gene knockout in human diploid cells.

The advance of CRISPR/Cas9 technology has enabled us easily to generate gene knockout cell lines by introducing insertion-deletion mutations (indels) at the target site via the error-prone non-homologous end joining repair system. Frameshift-promoting indels can disrupt gene functions by generation of a premature stop codon. However, there is growing evidence that targeted genes are not always knocked out by the indel-based gene disruption. Here, we established a pipeline of CRISPR-del, which induces a large chromosomal deletion by cutting two different target sites, to perform 'complete' gene knockout efficiently in human diploid cells. Quantitative analyses show that the frequency of gene deletion with this approach is much higher than that of conventional CRISPR-del methods. The lengths of the deleted genomic regions demonstrated in this study are longer than those of 95% of the human protein-coding genes. Furthermore, the pipeline enabled the generation of a model cell line having a bi-allelic cancer-associated chromosomal deletion. Overall, these data lead us to propose that the CRISPR-del pipeline is an efficient and practical approach for producing 'complete' gene knockout cell lines in human diploid cells.

Insertion/deletion (Indel) variant of the goat RORA gene is associated with growth traits.

RAR related orphan receptor A (RORA), which encodes the retinoid-acid-related orphan receptor alpha (RORα), is a clock gene found in skeletal muscle. Several studies have shown that RORα plays an important role in bone formation, suggesting that RORA gene may take part in the regulation of growth and development. The purpose of this research is to study the insertion/deletion (indel) variations of the RORA gene and investigate the relationship with the growth traits of Shaanbei white cashmere (SBWC) goats. Herein, the current study identified that the P4-11-bp and P11-28-bp deletion sites are polymorphic among 12 pairs of primers within the RORA gene in the SBWC goats (n = 641). Moreover, the P11-28-bp deletion locus was significantly related to the body height (p = 0.046), height at hip cross (p = 0.012), and body length (p = 0.003). Both of P4-11-bp and P11-28-bp indels showed the moderate genetic diversity (0.25

ARPIP: Ancestral Sequence Reconstruction with Insertions and Deletions under the Poisson Indel Process.

Modern phylogenetic methods allow inference of ancestral molecular sequences given an alignment and phylogeny relating present-day sequences. This provides insight into the evolutionary history of molecules, helping to understand gene function and to study biological processes such as adaptation and convergent evolution across a variety of applications. Here, we propose a dynamic programming algorithm for fast joint likelihood-based reconstruction of ancestral sequences under the Poisson Indel Process (PIP). Unlike previous approaches, our method, named ARPIP, enables the reconstruction with insertions and deletions based on an explicit indel model. Consequently, inferred indel events have an explicit biological interpretation. Likelihood computation is achieved in linear time with respect to the number of sequences. Our method consists of two steps, namely finding the most probable indel points and reconstructing ancestral sequences. First, we find the most likely indel points and prune the phylogeny to reflect the insertion and deletion events per site. Second, we infer the ancestral states on the pruned subtree in a manner similar to FastML. We applied ARPIP (Ancestral Reconstruction under PIP) on simulated data sets and on real data from the Betacoronavirus genus. ARPIP reconstructs both the indel events and substitutions with a high degree of accuracy. Our method fares well when compared to established state-of-the-art methods such as FastML and PAML. Moreover, the method can be extended to explore both optimal and suboptimal reconstructions, include rate heterogeneity through time and more. We believe it will expand the range of novel applications of ancestral sequence reconstruction. [Ancestral sequences; dynamic programming; evolutionary stochastic process; indel; joint ancestral sequence reconstruction; maximum likelihood; Poisson Indel Process; phylogeny; SARS-CoV.].

Reversal and Indel Distance With Intergenic Region Information.

Recent works on genome rearrangements have shown that incorporating intergenic region information along with gene order in models provides better estimations for the rearrangement distance than using gene order alone. The reversal distance is one of the main problems in genome rearrangements. It has a polynomial time algorithm when only gene order is used to model genomes, assuming that repeated genes do not exist and that gene orientation is known, even when the genomes have distinct gene sets. The reversal distance is NP-hard and has a 2-approximation algorithm when incorporating intergenic regions. However, the problem has only been studied assuming genomes with the same set of genes. In this work, we consider the variation that incorporates intergenic regions and that allows genomes to have distinct sets of genes, a scenario that leads us to include indels operations (insertions and deletions). We present a 2.5-approximation algorithm using the labeled intergenic breakpoint graph, which is based on the well-known breakpoint graph structure. We also present an experimental analysis of the proposed algorithm using simulated data, which showed that the practical approximation factor is considerably less than 2.5. Furthermore, we used the algorithm in real genomes to construct a phylogenetic tree.

Goat CLSTN2 gene: tissue expression profile, genetic variation, and its associations with litter size.

Calsyntenin-2 (CLSTN2) is involved in cell proliferation, differentiation, cell death, tumorigenesis, and follicular expression. Although CLSTN2 has been identified as a potential candidate gene for sheep prolificacy, no studies have been done on its effect on goat prolificacy. The purpose of this study was to identify mRNA expression and genetic variation within goat CLSTN2, and its association with prolificacy. Herein, we uncovered significant differences in mRNA levels of the CLSTN2 gene in different tissues in female goats (p < 0.01), including ovary tissue. Nine putative indels were designed to investigate their correlation to litter size, but only one 16-bp deletion was discovered in female Shaanbei white cashmere goats (n = 902). We discovered that a 16-bp deletion within the CLSTN2 gene was significantly correlated with first-born litter size (p = 0.0001). As shown by the chi-squared test, the genotypic II of single-lambs and multi-lambs was dramatically higher than with genotype ID (p = 0.005). Our findings suggest that indel within the CLSTN2 gene is a candidate gene affecting prolificacy in goats and may be used for Marker Assisted Selection (MAS) in goats.

Engineering indel and substitution variants of diverse and ancient enzymes using Graphical Representation of Ancestral Sequence Predictions (GRASP).

Ancestral sequence reconstruction is a technique that is gaining widespread use in molecular evolution studies and protein engineering. Accurate reconstruction requires the ability to handle appropriately large numbers of sequences, as well as insertion and deletion (indel) events, but available approaches exhibit limitations. To address these limitations, we developed Graphical Representation of Ancestral Sequence Predictions (GRASP), which efficiently implements maximum likelihood methods to enable the inference of ancestors of families with more than 10,000 members. GRASP implements partial order graphs (POGs) to represent and infer insertion and deletion events across ancestors, enabling the identification of building blocks for protein engineering. To validate the capacity to engineer novel proteins from realistic data, we predicted ancestor sequences across three distinct enzyme families: glucose-methanol-choline (GMC) oxidoreductases, cytochromes P450, and dihydroxy/sugar acid dehydratases (DHAD). All tested ancestors demonstrated enzymatic activity. Our study demonstrates the ability of GRASP (1) to support large data sets over 10,000 sequences and (2) to employ insertions and deletions to identify building blocks for engineering biologically active ancestors, by exploring variation over evolutionary time.

Discovery of DNA polymorphisms via genome-resequencing and development of molecular markers between two barley cultivars.

KEY MESSAGE: Genome resequencing uncovers genome-wide DNA polymorphisms that are useful for the development of high-density InDel markers between two barley cultivars. Discovering genomic variations and developing genetic markers are crucial for genetics studies and molecular breeding in cereal crops. Although InDels (insertions and deletions) have become popular because of their abundance and ease of detection, discovery of genome-wide DNA polymorphisms and development of InDel markers in barley have lagged behind other cereal crops such as rice, maize and wheat. In this study, we re-sequenced two barley cultivars, Golden Promise (GP, a classic British spring barley variety) and Hua30 (a Chinese spring barley variety), and mapped clean reads to the reference Morex genome, and identified in total 13,933,145 single nucleotide polymorphisms (SNPs) and 1,240,456 InDels for GP with Morex, 11,297,100 SNPs and 781,687 InDels for Hua30 with Morex, and 13,742,399 SNPs and 1,191,597 InDels for GP with Hua30. We further characterized distinct types, chromosomal distribution patterns, genome location, functional effect, and other features of these DNA polymorphisms. Additionally, we revealed the functional relevance of these identified SNPs/InDels regarding different flowering times between Hua30 and GP within 17 flowering time genes. Furthermore, we developed a series of InDel markers and validated them experimentally in 43 barley core accessions, respectively. Finally, we rebuilt population structure and phylogenetic tree of these 43 barley core accessions. Collectively, all of these genetic resources will facilitate not only the basic research but also applied research in barley.

Single-cell genome sequencing of human neurons identifies somatic point mutation and indel enrichment in regulatory elements.

Accurate somatic mutation detection from single-cell DNA sequencing is challenging due to amplification-related artifacts. To reduce this artifact burden, an improved amplification technique, primary template-directed amplification (PTA), was recently introduced. We analyzed whole-genome sequencing data from 52 PTA-amplified single neurons using SCAN2, a new genotyper we developed to leverage mutation signatures and allele balance in identifying somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) in PTA data. Our analysis confirms an increase in nonclonal somatic mutation in single neurons with age, but revises the estimated rate of this accumulation to 16 SNVs per year. We also identify artifacts in other amplification methods. Most importantly, we show that somatic indels increase by at least three per year per neuron and are enriched in functional regions of the genome such as enhancers and promoters. Our data suggest that indels in gene-regulatory elements have a considerable effect on genome integrity in human neurons.

Insertions and deletions as phylogenetic signal in an alignment-free context.

Most methods for phylogenetic tree reconstruction are based on sequence alignments; they infer phylogenies from substitutions that may have occurred at the aligned sequence positions. Gaps in alignments are usually not employed as phylogenetic signal. In this paper, we explore an alignment-free approach that uses insertions and deletions (indels) as an additional source of information for phylogeny inference. For a set of four or more input sequences, we generate so-called quartet blocks of four putative homologous segments each. For pairs of such quartet blocks involving the same four sequences, we compare the distances between the two blocks in these sequences, to obtain hints about indels that may have happened between the blocks since the respective four sequences have evolved from their last common ancestor. A prototype implementation that we call Gap-SpaM is presented to infer phylogenetic trees from these data, using a quartet-tree approach or, alternatively, under the maximum-parsimony paradigm. This approach should not be regarded as an alternative to established methods, but rather as a complementary source of phylogenetic information. Interestingly, however, our software is able to produce phylogenetic trees from putative indels alone that are comparable to trees obtained with existing alignment-free methods.

Development and validation of a forensic multiplex InDel assay: The AGCU InDel 60 kit.

Marker sets based on insertion/deletion polymorphisms (InDels) combine the characteristics of both short tandem repeats (STRs) and single nucleotide polymorphisms and have served as effective complementary or stand-alone systems for human identification in forensics. We developed a novel multiplex amplification detection system, designated the AGCU InDel 60 kit, containing 57 autosomal InDels, 2 Y-chromosomal InDels, and the amelogenin locus and validated the kit in a series of studies, which included tests of the PCR conditions; tests for sensitivity, species specificity, reproducibility, stability, and mock case samples; degradation studies; and a population study. The results indicated that the AGCU InDel 60 kit was accurate, specific, reproducible, stable, and robust. Complete DNA profiles were obtained even with 125 pg of human DNA. In tests of artificially degraded samples, we found that the number of alleles detected by the validated kit was considerably greater than that detected by the STR-based AGCU 21+1 kit, even as the degree of degradation increased. Additionally, 564 unrelated individuals from three Han groups were investigated using this novel system, and the values of combined power of discrimination and combined power of exclusion were not less than 1-4.9026 × 10-24 and 1-3.1123 × 10-5 , respectively. Thus, the results indicated that the novel kit was more powerful than the previous version of the InDel kit (the AGCU InDel 50 kit). Our results suggest that the AGCU InDel 60 kit can serve as an efficient tool for human forensics and a supplementary kit for population genetics research.

Learning the statistics and landscape of somatic mutation-induced insertions and deletions in antibodies.

Affinity maturation is crucial for improving the binding affinity of antibodies to antigens. This process is mainly driven by point substitutions caused by somatic hypermutations of the immunoglobulin gene. It also includes deletions and insertions of genomic material known as indels. While the landscape of point substitutions has been extensively studied, a detailed statistical description of indels is still lacking. Here we present a probabilistic inference tool to learn the statistics of indels from repertoire sequencing data, which overcomes the pitfalls and biases of standard annotation methods. The model includes antibody-specific maturation ages to account for variable mutational loads in the repertoire. After validation on synthetic data, we applied our tool to a large dataset of human immunoglobulin heavy chains. The inferred model allows us to identify universal statistical features of indels in heavy chains. We report distinct insertion and deletion hotspots, and show that the distribution of lengths of indels follows a geometric distribution, which puts constraints on future mechanistic models of the hypermutation process.

Efficient multinucleotide deletions using deaminase-Cas9 fusions in human cells.

CRISPR/Cas9 system is a robust genome editing platform in biotechnology and medicine. However, it generally produces small insertions/deletions (indels, typically 1-3 bp) but rarely induces larger deletions in specific target sites. Here, we report a cytidine deaminase-Cas9 fusion-induced deletion system (C-DEL) and an adenine deaminase-Cas9 fusion-induced deletion system (A-DEL) by combining Cas9 with rat APOBEC1 (rA1) and TadA 8e, respectively. Both C-DEL and A-DEL improve the efficiency of deletions compared with the conventional Cas9 system in human cells. In addition, the C-DEL system generates a considerable fraction of predictable multinucleotide deletions from 5'-deaminated C bases to the Cas9-cleavage site and increases the proportion of larger deletions at the target loci. Taken together, the C-DEL and A-DEL systems provide a practical strategy for producing efficient multinucleotide deletions, expanding the CRISPR/Cas9 toolsets for gene modifications in human cells.